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ABSTRACT 



We have investigated the formation of bound state of two electrons in different kind 
of lattices using a t-J-U model. In the model hopping parameter t tries to delocalize 
the electrons where as pairing of electrons comes via Heisenberg exchange integral J 
and hence it becomes necessary to calculate the threshold value of J, viz. J c , for which 
formation of bound states between two electron system is possible. The analysis is 
repeated for one dimensional chains, two dimensional square lattices, two leg ladders 
and three dimensional simple cubic lattices. Further we calculated the bound state 
energies for J > J c . Also we have tried to shed some light to the symmetry operation 
of the lattices to understand the characteristic of two electron pairing. 
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1 INTRODUCTION: 



The discovery of superconductivity in 1911 by Kammerlingh Onnes had given a new di- 
mension in Condensed Matter Physics research. From the very beginning the scientific 
and commercial potential of superconductors had been well understood by the commu- 
nity. So as the days progressed different exciting features of superconductivity started 
to come to light and in 1987 with the discovery of high temperature superconductors 
the field of interest is further broadened. 

Normal superconductors, which are usually good metals, are quite well understood 
by BCS theory where electron electron interaction is mainly controlled by phonons. 
But in high temperature superconductors the interaction picture among the electrons 
are still not clear. 

Here our motivation is to study strongly correlated systems because it is well agreed 
that the origin of high temperature superconductivity is purely from electronic inter- 
action, as for example here we are interested to study two dimensional square lattice 
which is analogous to the C11O2 planes in a high temperature superconductor. 

There are several models to study many particle systems and we are here using 
t-J model for studying electron pairing in metals with t denoting kinetic energy and 
J denoting Heisenberg exchange integral. We are starting with the assumption that 
the parent compounds are quite well represented by the antiferromagnetic Heisenberg 
model with localized electrons of spin 1/2 occupying a lattice point and coupled by an 
exchange integral J. Doping is assumed to remove electrons thereby producing "holes" 
which are mobile because neighboring electrons can hop to the hole site with amplitude 
t. It has been shown earlier that pQ for the t-J model, dilute holes in a antiferromagnet 
are unstable against phase separation into a hole rich and a no hole phase. It can argued 
that there exists a critical value of J c , such that when spin exchange integral J exceeds 
J c , the hole rich phase has no electrons and for J slightly less than J c the hole rich 
phase is a low density superfluid of electron pairs. 

Here we are investigating the critical value of J for which pairing of electron is 
possible in different kind of lattices. On the due course in the earlier semester we had 
concentrated on one dimensional chains. One dimensional analysis always carries a 
significant importance for its relative simplicity, also if one dimensional lattice features 
are well understood then it becomes relatively easy to understand the higher dimen- 
sional lattices. Then we had taken one step ahead to two dimensional square lattice. 
As we have stated the features of the two dimensional lattice has become very im- 
portant after the discovery of high temperature superconductors in ceramic materials. 
The structural feature of ceramic materials is Cu02 plane which are the main source 
for high temperature superconductivity. So it is important to understand the electron 
correlation in two dimensional lattices. Further we had extended our study on two 
leg ladder lattice. Since for the one-dimensional chain system quite a few things are 
known exactly. One approach to tackle the superconducting cuprates is to investigate 
the quasi one-dimensional lattices known as 'spin ladder' structures, which are strips 
of square lattice with a finite width and infinite length. An example of spin ladder 
lattice system is SV n _ 1 Cu n+1 2 n For these above mentioned lattices we had verified 



2 



the critical value of J in square lattice which is 2t[2j. Also we had investigated the 
critical value for one dimensional chain like lattice and two leg ladder where again we 
land up with the same result as J = It. 

In this semester we were curious to look at the critical value of J in a simple cubic 
lattice. To understand the electronic interaction of the heavy-fermionic superconduc- 
tors such as CeCuiSi?,, one needs to deal with different three dimensional crystal 
structures. So it is a natural to look at a three dimensional structure and try to find 
out the critical value of J for which pair formation among two electrons is possible. So 
we calculated J c for a simple cubic lattice and this constitutes a new and central result 
in this project work. 



2 MANY PARTICLE SYSTEM: 

2.1 Second Quantization: 

To understand the many body theory in condensed matter physics the essential tech- 
nique is the method of second quantization. Soon after the foundation of quantum 
theory, the formalism of creation and annihilation operator (second quantization) was 
introduced. The physics of creation and annihilation operators can be explained in a 
better way from relativistic quantum field theory. So second quantization is is nothing 
but a alternative formulation of quantum mechanics, the creation and annihilation 
operators are nothing but a tool that permits different process such as creation and 
annihilation of operators. Such process can not be discussed in the framework of 
Schrodinger equation[3]. 

We know that quantum mechanical wave function which represents a collection of elec- 
trons is antisymmetric with respect to the operation which exchanges the space and 
spin coordinates of any two electrons. Thus, if ?/>(l,2,3, ...,N) is an N electron wave 
function, and if Pij is the operator which exchanges the coordinates if electron i and 
electron j , then ^ = — ^ 

J. C. Slater has introduced a method to represent such many electron wave functions. 
We need to begin with an orthonormal set of one-electron functions: 1; 2 , where 
/ ^1^(1)^(1) = 5^3] , where 



A 



<Mn) <fa.(r N ) 



>N 



(n) 4>N(r N ) 



Such a function is called "Slater determinant." Since interchange of any two columns 
of a determinant changes its sign, A is antisymetric with respect to the exchange op- 
erator P^. 

We can also represent antisymmetric many electron wave function in a different man- 
ner. Let us now define a set of electron creation operators, &J,&2,...,, corresponding 
to the one electron spin orbitals, 0i, 02, When the creation operator &j acts on 
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an N-electron state, it produces an (N+l) electron state by creating an electron in 
the spin-orbital <f>j. We used to denote no electron state or "vacuum state" as |0 >. 
Similarly we define annihilation operator as b y 
The commutation relations of the operators as follows: 



6j-|0 >= 0, < 0\b] = 



\bj\N >=\N-1 >, b)\N >= 
bib] = Sij 



u i u 3 1 "3"i 



b\b) + b)b\ = 
hbj + bjbi = 



Now we can also represent Slater determinant with creation and annihilation 

0i(ri) <f>i(r 2 ) 

02 02(^2) 



operators, ^ 



Also ftM-lO >= -b]b]\0 > 



= 10102 



t^ti 



3% I 



Now we can write kinetic energy operator as 

f = Y.< i \T\j>b\b j 
Similarly, the potential energy operator can be written as, 

V = \ E <ij\V(r 1 ,r 2 )\kl>b\b]b l b k 

i,j,k,l 



Now if we define field operator as 



tt(r) = E^0™(r) 

i 

^W = E^0LW 



then 



< i\f\j >= J d?r<t> n ?{r)T{r)<t> ni {r) 

<ij\V\kl >= J J d 3 rid 3 r 2 0T. (n)0^. (r 2 )V\n - r 2 |0 ni (ri)0 ni (r 2 ) 

If we consider Bloch wave function i.e 0(r) = ^^Mfc n (r), where Uu n (j) signifies particle 
is in periodic potential the Hamiltonian will be, 

H = ^2,e k b\ a b k(T + - E Vqb\- q ^bki+ q ,obk'<T'bko 
this is second quantized Hamiltonian in many electron system. 

As for example if we consider e k = —t J2a e lk <5 , where t is hopping strength and 5 is 
nearest neighbour where a electron can hop. If our system is a two dimensional square 
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lattice with lattice parameter a then, 5 = ±xa, iya, 
So our energy dispersion relation will be then 

e k = -t(e ±ikx + e ±lky ) = -2t(cosk x + cosk y ) 

2.2 Theoretical Models in Many Particle Systems: 

2.2.1 Tight Binding Model: 

In tight binding model the Hamiltonian describes the kinetic energy (hopping) of elec- 
trons for nearest neighbour pairs. 

H = t (c^Cjc + Cj^Cia-) 

<i,j>,<j 

In this approximation we consider the wave function of the electrons are sharply lo- 
calized neglecting any overlap between them and they are confined in the lattice sites 
by an infinite potential barrier [1]. From here on we like to fix the notation of creation 
and annihilation operator as c$, cj respectively. 

2.2.2 Heisenberg Model: 

The simplest model in quantum many body theory is isotropic spin half Heisenberg 
chain. The Heisenberg Hamiltonian is given by: 

H = J £ (SSj - \) 

<i,j> 

where Sj is a local spin variable at j state. For antiferromagnets J > 0. The sum 
is over the distinct nearest neighbours. For spin 1/2 particle the spins are represented 
by Pauli's spin matrices. 

2.2.3 Hubbard Model: 

The Hubbard model describes the strongly correlated electron systems. The model in 
more than one dimension has not been solved. In describing the C11O2 planes in high 
temperature superconductivity Hubbard model is a good starting point [5]. The basic 
ingredients of Hubbard model are: 

• The kinetic energy (electron hopping) delocalizes the electron in Bloch state, 
leading to metallic behavior. [U] 

• The electron electron interaction (approximated by onsite Coulomb interaction) 
wants to localize the electron on to sites. [6] 
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The Hubbard model contains only one orbital per site and is defined as (considering 
nearest and next nearest neighbour interactions) [5J: 

h = -t ( c tcja + c) a c ia ) - 1' J2 (4a c i'<r + 4 a c la ) + uJ2 n n n n 

<i,j>,a <i,i'>,a i 

If we only take into account the nearest neighbour interaction it reduces to: 
H = —t (clcja + c) a c la ) + UJ2 ^T n a 

<i,j>,cr i 

where cj^ are creation operators and rii a = cl^c^ are occupation number operator. By 
means of on site Coulomb U the singlet band of the Hubbard model is split into a lower 
(LHB) and an upper Hubbard band (UHB). But the validity of the three band model 
to the single band model is still controversial It has been questioned whether the 
strong coupling version of the Hamiltonian, i.e the t- J model, is appropriate to describe 
correctly the low energy physics of the original three band model. In the large U limit 
and at half filling (one electron per site) the Hubbard ladder is equivalent to the spin 
1/2 Heisenberg ladder [7]. 

h — j Sj • Sj 

<i,j> 

When the Hubbard ladder is doped with holes away from half filling, its strong coupling 
description is modified from the Heisenberg model to the t-J model with the constrain 
of no doubly occupied sites [TJ. 

H = -t ]T {c\ a c ja + c) a c ia ) + J (Si • Sj - -rnuj) 

<i,j>,a <i,j> ^ 

2.2.4 t-J-U Model 

The most important local interactions in a doped antiferromagnet are well represented 
by the large U Hubbard model, the t-J model, and their various relatives. To be 
concrete we will focus on the t-J-U model. The t-J-U Hamiltonian is written as: 

H=-t £ £(cU CT + Cj W) + JE(Sr Sj - ^) + UJ2 nnrm 

<i,j> <* <i,j> i 

It is a close variant of the familiar t-J model where the sites of the lattice is strictly 
prohibited against double occupancy, i.e the doubly occupied sites are projected out. 
The "no double occupancy" restriction is suitably achieved by using "constrained" 
fermionic operators. The same physics can be achieved by using a t-J-U model given 
by the limit U — > oo we go back to simple t-J model. The exchange integral J arises 
through virtual processes where in the intermediate state has a doubly occupied site, 
producing an antiferromagnetic coupling. Dopping is assumed to remove electrons 
thereby producing a "hole" or missing spin which is mobile because neighbouring 
electrons can hop into its place with amplitude t. So among the other models to 
study the correlation effects in the high temperature superconductors this model is 
simplest one and gives a exact critical value of J(J C ) such that J > J c two electrons 
can form a two particle bound state. 
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2.3 Construction of Equation of Motion (EOM) for Two Par- 
ticle System: 

The system consisting of two electrons the wave function can be written as 

I*) = E $ (^2)cf lT ijO > (1) 

and the model Hamiltonian is 

H = -t £ XXcU, + 4c*) + J £ - ^) + C/£ n, T n a (2) 



where |0) denotes the vaccum state. For a two body problem the ground state isa 
singlet i.e &(i 1 ,i 2 ) = $(«2«i) and we know that H\^} = E\^) so the equation of 
motion can be written as 

H E E( c U.+sW)+j x: (^.^-^)+f/5:n lT n li ]E $ (^^ 2 )4 1 T4ilo>=^E $ (^^ 2 ) 

<i,j> cr <i,j> i U,i2 U,i2 

(3) 

*» = -* E Efe + iOE ^i^2)4 lT 4jo) 

= -*EE Hii,i2)(cl^c\^clj0)+c^ 

<i,j> h ,12 

= £ $0', i 2 )ct T cJ 2 jO> + - £ Sfa, i)44 lT |0) + 

*l*> = -^E^kWio) - $(ii,i)44 lT |o)] (4) 



<i,j> iiM 



U,»2 



« = iEE ^(<i,<2)(s, • sj - ^)c! lT c! 2i |o) 

<i,j> 11,12 

= j E $ (^) c utiil > ( 6 ) 

U,12 

£$(ii,i 2 ) = J2[thMj,i2)+t i2j $(i 1 ,j)} + [U5 lul2 - J il)i2 ]$(M2) (7) 
Fourier transform of the equation yields 

E$(k l5 k 2 ) = [t(k 1 )+t(k 2 )]$(k 1 , k 2 )+^ X! $ ( k i + k, k 2 - k)-j- £ J(k)$(k 1 - k, k 2 + k) 

(8) 



where, 

*(ki,k 2 ) = ^ E $ (ii,i 2 )e- ikl - r, i- ika - r '» 

*(k) = ^ E^' e ~ k ' (ri_rj) = -2t(coak x + cosfc y ) 

J(k) = 2J(cos k x + cos k y ) 

taking the lattice constant 1. Since the system is translationally invariant, the total 
momentum can be used to specify its eigenstates. let us define Q = k x + k 2 ,q = 
|(kx — k 2 ), and $(k!,k 2 ) = $Q(q)then we obtain 



®Q\<1) ~ ^ T7q— x T7q x I 9 ) 



U S fc $q(k)-^S fc J(q-k)$ Q (k) 
E-t(% +q)-t(f -q) 

This is the starting point of our analysis. 



3 CALCULATIONS: 

3.1 One Dimensional Chains: 



Figl 

For singlet pairing we can take Q=0 so from eq. (9) can be decoupled so that we can 
write, 

Co = UC I - 2JI X C X (10) 
C x = UCqI x — 2JI XX C X (11) 



where 



^o = ^E $ o(q) 



N , 



C * = T7E COS ^ $ o(k) 



N k 

j = 1e " 

N ^ E + 4tcosq x 

j 1 COSftc 



^ xx ~ N ^ 



cos 2 q 1 



E + At cos q x 



eq. (11) and (12) can be written in matrix form as follows: 

(UI -1 -2JI X \( C \ 

ul -2ji xx -l la 







;i2) 



In eq. (12) unique solution of C and C x will exit if and only if the determinant of the 
coefficient is zero. 

' UI -1 -23 L 



Now solving eq. (13) for J 



UL 



-23 



-23L 



1 







1-UI 



(13) 



(14) 



P x U + I xx (l-Uh) 
Using the lattice symmetry we can write 

j _ 1 1 
4 v / a ^I 

1 a + \J a 2 — 1 
* ~ 4 v^^T 

1 a(a + \Ja 2 — 1) 
J « = ~4 7^=1 

eq. (15), (16), (17) can be written is a more simple form as 

4 = — - — io (15) 

E . . 

hx = ( 16 ) 

substituting these values with the limit as I — > oo and [/ — > oo in eq. (14) we get 
Jc = 2t 



Bound state energy of electron in one dimension 

If we think in terms of the energy band in one dimension it is 8t for our system, so the 
bound state energy E of the two particle system can be obtained via equation (14), 
(15), (16) is we can write 

- = --- (17) 
3 tl t K J 

Taking t — 1 we can write the energy equation as 



r, 1 8t 

E = 

/o 3 



(18) 
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3.2 Two Dimensional Square Lattice: 




Fig 3 

In a similar fashion as we have done in the previous section we can write the following 
equations: 

C = UI C - 2JI X C X - 2JI y C y (19) 
C x = UCqI x — 2JC X I XX — 2JCyI xy (20) 

Cy = U C()Iy — 2JC X I X y ~ 2JCyIyy (2l) 



where 



Co = ^E$o(q) 

q 

c * = ^Xl cos ^ $ o(k) 

k 

I 1 v 1 

N Y E + 4t(cos q x + cos q y ) 

j J_y- COSftc 

x ~ N Y E + 4t(cos q x + cos g„) 
11 



hx - , r E 



cos 2 q 1 



N fr* E + 4t cos g a 



cosq y 



i =-y— 

y N ^ E + 4t(cos + cos q y ) 



1 



cos q x cos g y 



'2/2/ 



N q E + 4t(cos + cos q y ) 

1 cos 2 g,, 
iV ^ E + 4t(cos ^ + cos q y ) 



eq. (20), (21), (22) can be written in a matrix form as follows: 

-23 h 



( UI - 1 

UI X 



-2JJ, 
-2J7. 



1 



-2J/ y ^ 








—2JI xy 








2JIyy — 1 / 









(22) 



In eq. (23) unique solution of Co, C x and C y will exit if and only if the determinant of 
the coefficient matrix is zero. 



Uh 



UI X 
Uh 



-2JI X 
-2JI XX 
-2JL 



-2JI y 
-2 J I 



xy 



xy 

2t/ I yy 1 



For an isotropic square lattice symmetry permits us to write I y 
C y = C x hence eq. (23) becomes 



(23) 

11 - 1 anrl 
± xi ± yy J -xx cinu. 



/ ui - 
ui x 

V Uly 



-23 h 



-233, 



2 J i. <r- qr 1 2 J i. ' i 



xy 



-231 



xy 



-23 h 




= 



(24) 



and our modified determinant will be, 



Uh 



-23h 



UI X 
UL 



-23h 



-23I X 
-231 



xy 



-231 



xy 



-23h 



(25) 



Now our motivation is to take U — > oo limit to project out the possibility of double 
occupancy. The value of I , I x , I xx , I xy after integration are as follows: 



r 

-In 



2 na 

-2\ 1 <„, n\n/2 -2> 



1(1 -a)K(^) i(q-2)n(g,^) 



net 



7T 



1 (o: 2 -2a + 2)K(=£) 1 (a - 2)II(|, gj 1 qjg(-l) 



7TO; 



7T 



7T 



(26) 

(27) 
(28) 



12 



l (l-a)K(^) i (l a -g + i)fl:(-2) laE(- 



XV 2 7T 2 TTCt 4 7T 



These expression can be written in a simple form as: 



(29) 



E 

= Ft ~ I'» (31) 

/.= M«(§) (32) 

So the lattice integrals are turned out in terms of complete elliptic integral of first kind 
(K(k)), second kind (E(k)), third kind (II (v, k)) respectively. II(i/, fc) can be evaluated 
through K(k), E(k), F(q,(f>), E(q, 0), where F(q,(f)) and E(q,(f>) are the incomplete 
Elliptic integrals of first and second kind respectively [H]. The nature of i^(/c) and -E(fc) 
is shown as follows: 




Fig 4 

So from the physical nature of the elliptic integral of the first kind we can conclude that 
it is of logarithmic diverging nature so to tackle this problem let us consider a = —2 — 5 
and expand the elliptic integrals w.r.t 5 then take 5 — > limit. Since the determinant 
is zero thus the coefficient of the diverging InS term should be equal to zero. 
So-i^ + 18±-i4 + K = 

32 7T 7T 4 7T 8 7T 

The solution of the quadratic equation for J= 2t and 7.32t. 
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Bound state energy of electron in two dimension 

From eq. (24) we can write 



-23 



1-UI 



2UI* + (I xx + I xy )(l-UI ) 



(33) 



From eq. (31), (32), (33), and (34) the bound state energy E of the two electron system 
can be obtained via 



m _ 7i f e 
T-2jt(!)"7 



Taking t — 1 we rewrite the equation as: 



E 



16 



v 7r K (a) ' 



(34) 



(35) 



0.13 



0.11 



0.09 



Plot 2 

Bound State Energy vs. J 



LU 



0.07 



0.05 



0.03 



0.01 



3.3 



3.7 



3.9 



Fig 5 



14 



3.3 Two Leg Ladder: 



■ ^ — i 

( 






\ 

u 


u 




t 







Fig 6 

So a ladder like lattice involves a one dimensional wave vector integral (viz over 5 X ) 
rather that two dimensional integral (over q x and q y ). More precisely the lattice inte- 
grals appearing in the calculation of bound states are expressed as, 

e4 e ^ f_jo* 



^ q y =0,2n 2 71 " 



The various lattice integrals are as follows: 

dx 













' — 7T 




f'TV 








I—7Y 




f'TY 


I XX = ^ 






J — 7V 



■f 1 + COS ^ 

cos (fcrfx 



+ 



1 

4 7-7T ct - 1 



dx 



+ 1 + COS ^ 

cos 2 

+ 1 + cos q x 



+ 



1 Z" 71 " CO' 

4 7-7T a - 1 



+ cos q x 
cos (fccte 



+ 



+ cos q x 
cos 2 g x rfa; 



1 r n cos 
4 7-7T a - 1 



+ cos q x 



where a 

from eq. (13) with t — 1 and [/ — > oo we can write, 



if 



-13 U 



I x 2,3 1 xx 1 
expanding the determinant we obtain 



2 3I 2 X - 2 J/q/xx - Io = 



(36) 
(37) 
(38) 



(39) 



(40) 



We have substituted the integral values in eq. (42) and also a is replaced by = — 5 — 2. 
Then expanding the entire equation w.r.t 5 we got a term of 4=. This term will diverge 
as soon as we will take the limit as 5 — > 0. So collecting the coefficients of this term 
and equating them to zero the critical value of J is obtained which is 2t. 



Bound state energy of Two leg ladder 

Directly from the integrals one can find out the bound state energy for pairing of 
electrons with proper substitution. From the plot also we can verify that the minimum 
energy required for formation of bound state among two electrons in a ladder like 
lattice is 2, considering t=l. 
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Fig 7 
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3.4 Three Dimensional Cubic Lattice: 



u- 




Fig 8 

From eq.(9) considering that the lattice is isotropic (Q=0), we can write 

|E^o(k)4E^(q-kWk) 



$o(q) 



$o(q) 



u_ 

N 



E - 2f (q) 
E fe $o(k)4E^(q-kWk) 



E + 4t(cos q x + cos q y + cos q z ) 
This equation can be written as follows: 

Cq = U CqIq — 2 J C X I X — 2 J Cyly — 2 J C Z I Z 



where 



C x = UCqI x 2JC x I xx 1JCyI X y 2JC Z I X 

Cy = U C()Iy 2JC X Iy X 2JCyIyy 2JC Z Iy 

C z = U CqI z — 2JC X I ZX — 2JCyI zy — 2JC Z I Z 
Co = ^E $ o(q) 

q 

c * = ^H cos ^ $ o(k) 
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°y = T?E cos/c ?/ $ o(k) 



N 

tx, 

c * = ^E cos ^ $ o( k ) 

k 



1 



N q E + 4t(cos q x + cos q y + cos g z ) 



ix = — Y 

x AT ^ 



cosq x 



N q E + 4t(cos q x + cos g y + cos q z ) 



COS" ^ 



N g E + 4t(cos q x + cos q y + cos q z ) 



1 



cosq y 



N ^ E + 4t(cos + cos g y + cos q z ) 

1 cos g x cos g y 

^ N y E + 4t(cos g^ + cos q y + cos g z ) 



E 



cos^ ^ 



N q E + 4t(cos + cos g y + cos q z ) 

j = J_ cosg 2 

N g E + 4t(cos g^; + cos g y + cos q z ) 



cos 2 q z 



N g E + 4t(cos ga; + cos g y + cos g z ) 



cos q z cos g a 



N g E + 4t(cos ga; + cos q y + cos q z ) 

1 ^ cos g y cos q z 

N q E + 4t(cos g x + cos q y + cos g z ) 



eq. (44), (45), (46), (47) can be written in a matrix form as follows: 



/ UI - 1 -2 J/ w 



ui x 

Uly 



-2JJ. 
-2J/ 2 



-2J/ 2 

2JI Z x 

2t//yy 1 1JI yZ 



\ ( Co \ 

C x 



-2 J I, 



yz 



-2ji zz - 1 y 



v c, y 



= 



(47) 



J c /i can be found out by vanishing of the determinant by choosing an energy E to be 
slightly between noninteracting 2 electron band. 



UI - 1 -2JI X 



-2JI y 
-2JI. 



xy 



Uly 

UL 



-2 J I. 
-2JL 



xy 



-2JL 



yy 



-2JL 



yz 



-2JI Z 

2 J Izx 

2JI yz 

-2Jhz - 1 
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Since we are employing the constrain that there can not be any doubly occupied sites 
so considering U — > oo, our vanishing determinant be as follows: 



(4£ 



io 


-2JI X 


-2JI y 


-2JI Z 




2 J I n\ x 1 


2JI X y 




h 


2 J Ixy 


2t/ lyy 1 


2 J Iyz 


h 


2 J Izx 


2 J Iyz 


-2JI ZZ - 1 



Now because of lattice symmetry we can write 



1 yy 



J-xy ±yz 

Also these lattice integrals can be written as: 

1,1 E 



zz 
zx 



3 v 4t 
E 



At 



At uy 
I —21 

A x lL ' ± xy 



(49) 
(50) 



So we need to find out Iq and I xx or I xy . 



Solution of Lattice Integrals 

Previously mentioned lattice integrals are commonly known as lattice Green's function. 
The most general representation of lattice green's function is,[H] 



G(s, I, m, n) 



1 r n r n ft cos Ix cos my cos nzdxdydz 



7T° 



1 /*7T /*7T /*7T 

= G(s, 0, 0, 0) = - — 

7T d JO JO JO S — 



o jo jo s — (cos x + cos y + cos z) 

dxdydz 



(cos x + cos y + cos z) 

This integral defines a single valued analytic function G(s) in the complex s-plane cut 
along the real axis from -3 to +3. in most physical applications one usually requires 
the limiting behavior of the Green function G(s) as s approaches the real axis|lU]. 
Another way we can write as[TT] 



1 



G(s) 



1 1 

3^ 



o Jo 



dxdydz 



o s — i(cosx + cosy + cos 2) 



GAs) 



7T° JO JO 



dxdydz 



s — 3 (cos x + cos y + cos z) 



(51) 



The series representation of the Green's function is|llj: 

° n=0 
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(52) 



where j=l, 2, 3, 1 < \s\ < oo and 

Pn ] = A r r r[H^^es)] n de x de 2 de z 

TT JO JO Jo 

Now for simple cubic lattice Ai(6 l i, 6*2, 6*3) = |(cos#i + cos #2 + cos #3) The recurrence 
relation for simple cubic lattice coefficient p\ n is 

36(n + f) 3 pS + 2 " 2(2n, + f )(f0n 2 + lOn + 3)p£ + n(4n 2 - l)pg_ 2 = (53) 

from eq. (52) and (53) G\(s) can be written as a solution of the linear third order 
differential equation 

( s 2 _ l)(9s 2 _ + g s(9s 2 _ + 3(21s 2 _ 4) ^1 + QsGi = Q (54) 



So G\{s) can be written as, 



l-9{ 4 r 2^,„ 2 



°'M = ,(i-fl.(I + 3f) '; K( * l)1 (55) 



where, 



fee 3 



\ (i-0 3 (i + 30 



-1/2 / / — v 1/2 



and /«(A;i) is the complete elliptic integral of the first kind. 
Also 

Gi(l) = 3(f8 + f2 v / 2-f0v / 3-7v / 6)[-/t((2- ^(v 7 ^- v 7 ^))] 2 (56) 

7T 

So from this knowledge with proper adjustment of coefficients we can write: 

J ° = fk Gl(s) (57) 

Here we like to mention that elliptic integrals diverges for modulus of f , the physical 
significance of the divergence in our work is the electrons are on the energy band there- 
fore in our earlier analysis we had taken care of that divergence keeping our calculation 
just below the band, but here when we substitute our parameters as electrons are on 
the energy band we found the modulus is not f therefore there is no any sharp diver- 
gence of elliptic integral rather there is a local maxima, and we are returning with a 
numerical value of the integral on the band. Its a striking feature of the elliptic integral 
for simple cubic lattice. So from this knowledge we can evaluate Iq. But one of the 
integral among I xx or I X y are still to be evaluated. To evaluate our integrals we had 
taken help from the recurrence relation of the Green function. With the knowledge of 
the recurrence relation for fee lattice[f2j and triangular lattice[T3], we have constructed 
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the recurrence relation for simple cubic lattice. For nearest neighbour (6 nn for S.C) 
interaction only the recurrence is as follows: 

G(l+l,m,n)+G(l-l,m,n)+G(l,m+l,n)+G(l,m-l,n)+G(l,m,n+l)+G(l,m,n-l) 

= 25[ 5 m0 5 n o - 2sG(l, m, n) (58) 

where, 

1 r n f n f w cos (I + l)x cosmy cosnz 

G{1 + l,m,n) = — / / / — ■ ■ -dxdydz 

it A Jo Jo Jo s + (cos x + cos y + cos z) 

and so on, Let 1=1, m=0, n=0 then 

G(2, 0, 0) + (7(0, 0, 0) + G(1, 1, 0) + G(1, -1, 0) +G(1, 0, 1) +G(1, 0, -1) = -2sG(l, 0, 0) 

Due to the symetric structure of the lattice (7(1,1,0) = (7(1,-1,0) = (7(1,0,1) = 
(7(1, 0, -1) Hence, (7(2, 0, 0) + (7(0, 0, 0) + 4(7(1, 1, 0) + 2s(7(l, 0, 0) = 
Adjusting the coefficients for our calculation, 

where, I xx = ^((7(2, 0, 0) + G(0, 0, 0)), J = £(7(0, 0, 0), 4 = £(7(1, 0, 0) Solving eq. 
(50) and (59) I xy is determined. 

hy = -| 4 " ^0 (60) 

Let us now summaries the lattice integral values: 

In = (7i(s), S = — 

12* U ; At 
El 

Ixy = -g^4 - ^0 ( 61 ) 
T — T — 21 

Substituting eq. (60) in eq. (48) and solving for J, the values are obtained as 2, 2, 
7.88. 



4 RESULTS AND DISCUSSION: 

In our entire work we were involved in finding out the critical values of J ( J c ) for which 
bounstate formation of two electron system is possible, and our results obtained are 
as follows: 
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Lattice Type 


J c Values 


One Dimensional Chains 


2t 


Two Dimensional Square Lattice 


2t, 7.32t 


Two Leg Ladder 


2t 


Three Dimensional Cubic Lattice 


2t, 2t, 7.88t 



Also we have plotted J c /t Vs. \E bs \/t for one dimensional chain, two dimensional 
square lattice and two leg ladder. 

For one dimensional chain like system has been already solved for J=2tp3, [15]. 
For square lattice the critical value of J=2t can seen in papers of Lin [HE]. In our 
work we are proposing that for two leg ladder lattice the critical value of J is 2t, and 
it is quite desirable because two leg ladder is a simplified form of n leg ladder system 
which is nothing but a 2 dimensional system. 

We are also proposing for the first time that for three dimensional lattice the critical 
value of J is 4. 

Let us now try to understand the significance of the multiple value of two dimensional 
and three dimensional lattice for the basis of group theory. 
Square Lattice 

We know that the symmetries of a square lattice is represented by C^v group. The 
character table for C^y group is as follows: 



Characters 


E 


c 4 


cl 


cl 


m x 


m y 








1 


1 


1 


1 


1 


1 


1 


1 


x (2) 


1 


-1 


1 


-1 


-1 


-1 


1 


1 


x m 


1 


-1 


1 


-1 


1 


1 


-1 


-1 




1 


1 


1 


1 


-1 


-1 


-1 


-1 




2 





-2 


















From the character table we can find out the irreducible representation of the group 
and the corresponding basis functions[16j. 



Irreducible Representation 


Basis Functions 


r(i) 


1 


r( 2 ) 


xy 


r(a) 


x 2 — y 2 


r( 4 ) 


xy(x 2 — y 2 ) 


r(s) 


(x,y) 



The possible pairing symmetries correspond to these irreducible representation and to 
the basis functions. In case of square crystals the singlet orders are called s, d x 2_ y 2, d xy 
and g. The corresponding basis functions are 1, x 2 —y 2 , xy, xy(x 2 y 2 ). It is reasonable 
to classify order parameters as "s-wave like" and " d- wave- like" . So in our analysis 
J=2t corresponds to s-wave and J=7.32t corresponds to d-wave[7]. 
Simple Cubic Lattice 

It is known that the simple cubic lattice corresponds to Oh group. The character table 
for Oh group is as follows [T7j: 
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Characters 


E 


8C 3 


6C 2 


6C 4 


3C 2 (= ci) 


i 


6S 4 


8S e 




6o- d 


Mg 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Mg 


1 


1 


-1 


-1 


1 


1 


-1 


1 


1 


-1 


Eg 


2 


-1 








2 


2 





-1 


2 





Tig 


3 





-1 


1 


-1 


3 


1 





-1 


-1 


T2g 


3 





1 


-1 


-1 


3 


-1 





-1 


1 


A lu 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


A 2u 


1 


1 


-1 


-1 


1 


-1 


1 


-1 


-1 


1 


E u 


2 


-1 








2 


-2 





1 


-2 





T\ u 


3 





-1 


1 


-1 


-3 


-1 





1 


1 


T2g 


3 





1 


-1 


-1 


-3 


1 





1 


-1 



So the irreducible representation and basis functions for singlet pairing are [7] 



Irreducible Representation 


Basis Functions 


ri(A lfl ) 


1 


r 2 (A 2g ) 


{k\ — ky)(ky — kD(kl — kD 


UEg) 


ow un un un un 

n = 2; 4 


r 4 (T l9 ) 


U U (un _ un\ u u (un _ un\ u U (Un _ un\ 

n = 2;4;6 


r 5 (T 29 ) 


U U un u u un u U un 
r ^x ru y ru z ' ^y^z^x 5 nj Z rh X rh y 

n = 0;2;4 



So in a cubic lattice, Ti(Ai g ) represents s-wave pairing, T 3 (E g ) (for n=2) and r 5 (T 29 ) 
(for n=0) represents d-wave pairing|7J. Also T 3 (E g ) is doubly degenerate. So from 
the evaluated value of of J, we can conclude that J=2t value corresponds to T 3 (E g ) 
representation and J=7.88t corresponds to T 1 (A 2g ) representation. 

5 CONCLUSION AND FUTURE PROSPECTS: 

The analytical result of J for simple cubic lattice is a new result. The result can be 
cross examined by bilayer consideration of the three dimensional cubic lattice. The 
bilayer consideration is analogous with the two leg ladder consideration. Physically 
it can be vilualised that, n number of layers togather is forming a three dimensional 
structure provided the separation between the layers is sufficiently small. So starting 
the calculation for bilayer one can extend the calculation to n number of layer which 
can give the similar result as of the three dimension. A natural extrapolation can be 
investigating superconductivity in these lattices and studying the phase diagram. 
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